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ABSTRACT. The thickness vibrations of a finitely deformed infinite periodic laminate made 
out of two layers of dielectric elastomers is studied. The laminate is pre-stretched by inducing 
a bias electric field perpendicular the the layers. Incremental time-harmonic fields superim- 
posed on the initial finite deformation are considered next. Utilizing the Bloch-Floquet theorem 
along with the transfer matrix method we determine the dispersion relation which relates the 
incremental fields frequency and the phase velocity. 

Ranges of frequencies at which waves cannot propagate are identified whenever the Bloch- 
parameter is complex. These band-gaps depend on the phases properties, their volume fraction, 
and most importantly on the electric bias field. Our analysis reveals how these band-gaps can 
be shifted and their width can be modified by changing the bias electric field. This implies 
that by controlling the electrostatic bias field desired frequencies can be filtered out. Represen- 
tative examples of laminates with different combinations of commercially available dielectric 
elastomers are examined. 

Keywords: dielectric elastomers; wave propagation; finite deformations; thickness vibra- 
tions; non-linear electroelasticity, band-gap, composite, Bloch-Floquet analysis. 

1. Introduction 

Band-gaps corresponding to ranges of frequencies at which waves cannot propagate earned 
the interest of the scientific community (e.g., Ziegler, 1977; Wang and Auld, 1985; Kushwaha 
et al., 1994; Gei et al., 2004, 2009, to name a few). These band-gaps (BGs) or stop-bands 
can be utilized for various applications such as filtering elastic waves, enhancing performance 
of ultrasonic transducers, supplying a vibrationless ambient when needed for industrial and 
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biomedical applications, and more. Several techniques are available for investigating wave 
propagation in composites, (e.g., for the plane-wave method in periodic composites Kushwaha 
et al., 1993; Tanaka and Tamura, 1998; Sigalas and Economou, 1996) and the effective medium 
methods (e.g., Sabina and Willis, 1988, for composites with random micro structre). As in this 
work we are interested in the appearance of BGs in infinite periodic laminates, we find it 
convenient to use the transfer matrix method (Wang and Auld, 1986), along with the Bloch- 
Floquet theorem (Kohn et al., 1972). Specifically, we consider waves superposed on a pre- 
deformed state due to bias electric field acting on an infinite periodic laminate with a repeating 
sequence of of dielectric elastomers (DEs) layers. The motivation for the choice of DEs as the 
composite constituents stems from their ability to undergo large deformations (Pelrine et al., 
2000), and to change their mechanical and dielectric properties when subjected to electrostatic 
fields (Kofod, 2008), thus effecting the way in which electroelastic waves propagate in the 
matter (Shmuel et al., 2011; Gei et al., 2011). As will be shown in the sequel, by properly 
tuning the bias electrostatic field different ranges of frequencies can be filtered out. 

The structure of this paper is as follows. Following the works of Dorfmann and Ogden 
(2005); deBotton et al. (2007); Dorfmann and Ogden (2010); Bertoldi and Gei (201 1); Rudykh 
and deBotton (2011); Ponte Castaneda and Siboni (2011) and Tian et al. (2012), the back- 
ground required for describing the static and dynamic responses of elastic dielectric laminates 
is revisited in section 2. Particularly, the governing equations for small fields superimposed 
on large deformations of electroactive elastomers are summarized. In this section we also 
introduce an adequate extension of the transfer matrix method and the Bloch-Floquet theo- 
rem required for treating the propagation of incremental electroelastic waves superimposed on 
finitely deformed DE laminates. Section 3 deals with the response an infinite periodic lam- 
inate with alternating layers, whose behaviors are governed by the incompressible dielectric 
neo-Hookean model (DH), to a fixed electric displacement field normal to the layers plane. 
Subsequently, the response to small harmonic perturbations superimposed on the finitely de- 
formed laminate is analyzed. The end result is given in terms of a dispersion relation relating 
the incremental fields frequency and phase velocity to the Bloch-parameter. Several examples 
are considered in section 4, where the influence of the morphology of the laminate, the contrast 
between the phases properties, and the bias electric field on the dispersion relation is examined. 
Finally, realizable tunable isolators are studied by considering combinations of commercially 
available DEs such as VHB 4910, fluorosilicone 730 and ELASTOSIL RT-625. 
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2. THEORETICAL BACKGROUND 

Consider a heterogeneous body occupying a volume region £Iq C M 3 made out of N per- 
fectly bonded different homogeneous phases. The external boundary of the body dQ.Q sepa- 
rates it from the surrounding space IR 3 \Qo> assumed to be vacuum. Each phase occupies a 
volume region (r = l,2,...,N) and its boundary is dCl^. Let # : £2o x y — > M 3 describe 
a continuous and twice-differentiable mapping of a material point X from the reference con- 
figuration of the body to its current configuration Q. with boundary dQ. by x = % (X,?). The 
corresponding velocity and acceleration are denoted, respectively, by v = % t and a = % tv and 
the deformation gradient is F = V x #. Vectors in the neighborhood of X are transformed into 
vectors in the current configuration via dx = FdX. The volume change of a referential volume 
element dV is given by dv = JdV, where dv is the corresponding volume element in the current 
configuration, and J = det (F) > due to material impenetrability. The conservation of mass 
implies that pi = Jp, where pL and p are the material mass densities in the reference and the 
deformed configurations, respectively. An area element dA with the unit normal N in the refer- 
ence configuration is transformed to a deformed area da with the unit normal n in the current 
configuration according to Nanson's formula NdA = jF T nda. The right and left Cauchy-Green 
strain tensors are C = F T F and b = FF T . 

The electric field in the current configuration, denoted by e, is given in terms of a gradient of 
a scalar field, termed the electrostatic potential. In free space an induced electric displacement 
field d is related to the electric field via the vacuum permittivity Gq, such that d = £oe. In 
dielectric bodies the connection between the fields is specified by an adequate constitutive 
relation. 

In terms of a 'total' stress tensor o the equations of motion read 

V-<r = pa, (1) 

where in order to satisfy the balance of angular momentum symmetry of o is required. The 
terminology 'total' stress is used to indicate that o accounts for both mechanical and electrical 
contributions. Thus, the traction t on a deformed area element with a unit normal n is given by 

an = t. (2) 

When ideal dielectrics are considered no free body charge is present, and Gauss's law takes the 
form 

Vd = 0. (3) 

Faraday's law is 

V x e = 0, (4) 
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when a quasi-electrostatic approximation is considered. A full description of the system should 
consider interactions with fields outside the body, henceforth denoted by a star superscript. 
Specifically, these are 



e e 



I 



(5) 
(6) 



where I is the identity tensor. In vacuum the governing Eqs. (3)-(4) for d* and e* are equiva- 
lent to Laplace's equation for the electrostatic potential. As a consequence, o* known as the 
Maxwell stress, is divergence-free. Across the outer boundary d£l, the electric jump conditions 
are 



(d-d*)-n=-w e , (e-e*)xn = 0, (7) 

where w e is the surface charge density. In order to formulate the jump in the stress we postulate 
a separation of the traction into a sum of a prescribed mechanical traction t m , and an electric 
traction t e induced by the external stress such that t e = <7*n. Hence, the stress boundary con- 
dition is 

(CT-CT*)n = t m . (8) 

The jump conditions across a charge-free internal boundary between two adjacent phases m 
and / are 

[[a]]n = 0, [[d]]-n = 0, [[e]]xn = 0, (9) 

where [[•]] = - denotes the jump of fields between the two phases. 

A Lagrangian formulation of the governing equations and jump conditions is feasible by a 
pull-back of the Eulerian fields 

P = /aF T , D = JF 1 d, E = F T e, (10) 

for the 'total' first Piola-Kirchhoff stress, Lagrangian electric displacement and electric field, 
respectively. Thus, the governing Eqs. (1), (3) and (4) transform, respectively, to 

V x -P = p L a, V x -D = 0, V x xE = 0. (11) 

The corresponding referential jump conditions across the external boundary 3Q.Q are 

(P-P*)N = t M , (D-D*)-N = -w £j (E-E*)xN = 0, (12) 

where t M &A = t m da and wsdA = w e da. Across the referential interfaces Eqs. (9) become 

[[P]]N = 0, [[D]]-N = 0, [[E]]xN = 0. (13) 
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Following Dorfmann and Ogden (2005), the constitutive relation is expressed in terms of an 
augmented energy density function (AEDF) *F with the independent variables F and D, such 
that the total first Piola-Kirchhoff stress and the Lagrangian electric field are derived via 

P =5F- e =5d (14) 

The first of Eq. (14) should be modified when considering incompressible materials, as the 
kinematic constraint yields an additional workless part of stress. The latter is accounted by in- 
troducing a Lagrange multiplier po which can determined only from the equilibrium equations 
together with the boundary conditions. Thus, the total first Piola-Kirchhoff is 

p= ot-" oF <15) 

Based on the framework developed in Dorfmann and Ogden (2010), we superimpose an 
infinitesimal time-dependent elastic and electric displacements x = x(K,t) and D(X,f), on 
the pre-deformed configuration. Herein, and in the sequel, we use a superposed dot to denote 
incremental quantities. Let the Eulerian quantities E,d, and e denote the push-forwards of 
increments in the first Piola-Kirchhoff stress, the Lagrangian electric displacement and electric 
fields, respectively, namely 

E=-^PF r , d=-FD, e = F r E. (16) 
In terms of these variables, the incremental governing equations read 

V-E = px», V-d = 0, Vxe = 0. (17) 

For an incompressible material the incremental constraint is 

V-x = tr(h)=0, (18) 

where h = Vx is the displacement gradient. Linearization of the incompressible material con- 
stitutive equations in the increments yields 

E = tfh + p h T -pol + M, (19) 
e = ^ T h + ^d, (20) 

where (^ T h) fc = ^y^y- The quantities and are the push-forwards of the referential 

reciprocal dielectric tensor, electroelastic coupling tensor, and elasticity tensor, respectively. In 
components form 

^ij = J F oU J ^0ap F p~j , 9§ijjk = Fja&oiapFpk > %jU = j F ja%iakp F ip , (21) 

where 

d 2x ¥ d 2x ¥ d 2x ¥ 

= dDjD? ^ = dF~dDp %iakP = dF^dF^- (22) 
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The incremental outer fields are given by 

d* = £ e\ (23) 
6* = e [e* <g> e* + e* ® e* - (e* • e*) I] . (24) 

Here again d* and e* are to satisfy Eqs. (3)-(4), where subsequently V ■ 6* = 0. 

The push-forward of the increments of Eqs. (12) gives the following incremental jump con- 
ditions across the external boundary in the current configuration 

[E-(i* + a*h T -(V-x)(T*]n = t m , (25) 
[d-d'*-(V-x)d* + hd*] n = -w e , (26) 
[e - e* - h V] x n = 0, (27) 

where t m da = ij^dA and w e da = w E dA. Similarly, the push forward of the increments of 
Eqs. (13) results in the following jumps across the internal boundaries 

[[E]]n = 0, [[d]]-n = [[e]]xn = 0. (28) 

Several techniques are available in the literature for tackling the problem of wave propaga- 
tion in periodic composites. In layered structures it is convenient to utilize the transfer-matrix 
method in conjunction with the Bloch-Floquet theorem (Adler, 1990). Herein we formulate an 
adequate adjustment for treating the propagation of incremental electroelastic waves superim- 
posed on a finitely deformed multiphase DE laminates. 

Consider a finitely deformed infinite DE-laminate made out of m-phase periodic unit-cells 
subjected to incremental perturbations harmonic in time and normal to the layers plane. A 
schematic illustration of the laminate is shown in Fig. 1. Let s n p \ p = l,..,m, denote an 
incremental state vector in the p-phase of the n th cell consisting quantities which are continuous 
across the interface between the two adjacent layers. In view of the linearity of the incremental 
problem it can be shown that the state vector at the top of the p-phase s|fl and the state vector 
at the bottom of the p-phase s^K (see Fig. 1) are related via a unique non- singular transfer 
matrix such that 

<») 

In terms of the state vectors, the continuity conditions across the interface are 

^\ = i p A l) - (30) 

tin) b{n) v ' 

(m) (1) 

Utilizing Eqs. (29)-(30) recursively the relation between s ; and s) is 

s (b) = j(b) s (b) _ j(b) (m-l) _ j(b)j(m-l) (m-l) _ _ j 
t{n) b{n) t(n) b(n) "' b ' ^ ' 

where T = n»=i ^ s m e total transfer matrix. 
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Figure 1. An illustration of an infinite laminate composed of periodic m- 
layers unit-cells. 

At one hand the continuity conditions across the interface between two successive cells n 
and n + 1 are 



„(l) _ M 



(32) 



At the other, the Bloch-Floquet's theorem states that the state-vectors of the same phase in 
adjacent cells are identical up to a phase shift in the form 



Jj>) _ „-ik B hJj>) 



(33) 



where < kg < n/h is the Bloch-Floquet wavenumber in first irreducible Brillouin zone, rep- 
resenting the smallest region where wave propagation is unique (Kittel, 2005). Substitution of 
Eq. (31) in (32), followed by employing the latter with Eq. (33), yields the eigenvalue problem 



det 



T _ e -ik B h l 



0, 



(34) 



where I is an identity matrix with the dimension of the total transfer matrix. The solution of 
Eq. (34) provides the dispersion relation describing the manner by which waves propagate in 
the laminate. 
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3. Thickness vibrations of an infinite periodic two-phase dielectric 

laminate 

Henceforth we restrict our attention to isotropic phases. Consequently, the phases AEDFs 
can be written in terms of the six invariants (Dorfmann and Ogden, 2005) 

/ 1= tr(C) = C:I, I 2 = j(tf-C:C), I 3 = det(C) = J 2 , (35) 
7 4e = D D, 7 5e = D CD, 7 6e = D C 2 D. (36) 

Specifically, we consider phases which are characterized by incompressible dielectric neo- 
Hookean (DH) model 

fa,/*) = ^- (/i -3) + ^I 5e , (37) 

where fi^ and = £q£^ denote the phases shear moduli and dielectric constants, respec- 
tively, and Er being the phases relative dielectric constant. The corresponding total stress 
is 

G (p) = ^p) h + J-a®d-p i ( j\ (38) 

and the current displacement and electric field are related via the isotropic linear relation 

d = e (p) e. (39) 

Calculation of the associated constitutive tensors yields, in components 

form, 

where 8jj are the components of the Kronecker delta. 

Consider a Cartesian coordinate system with unit vectors 11,12 and i3 along x\,X2 and X3, 
respectively. Let a DE laminate be composed of periodic two-phase unit-cells with thickness 
H along the the X2 direction, infinite along the x\ and X3 directions as illustrated in Fig. 2. 
For convenience we denote the phases a and b. The sum of the phases thicknesses equals the 
unit-cell thickness H = +H^ b \ such that = /H and = H W /H are the volume 
fractions of phase a and phase b, respectively. We assume that the composite is allowed to 
expand freely in the X{ and X3 directions. Of interest is the composite response when subjected 
to an overall electric displacement field d = d^i, where J2 = v (a 

In each phase we assume displacement fields compatible with the homogeneous uni-modular 
(on account of incompressibility) diagonal deformation gradients 



[F] (p) =diag 



j (p) j (p) fi (p)\ 1 
1 ' 2 ' I 1 2 



(41) 
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Figure 2. An illustration of an infinite periodic laminate, composed of alter- 
nating a and b phases, with initial thicknesses H^ a ' and H^ b \ respectively. The 
electric bias field is along the X2 direction. 

where henceforth p = a,b. In virtue of the perfect bonding between the phases, the stretch 
ratios in the x\ and xj directions in the phases are the same, i.e., 

-1 / f„i -1 



x[ a ' 



5 (a) 5 (a) 
Aj A 2 



(42) 



Eqs. (42) imply that A^^ = Aj , and henceforth we define A = A^ a . From the symmetry of 



(p) 

the problem in the x\ — xt, it follows that A 2 



X 2 . The corresponding stresses in each phase 



are 



°11 — °33 



{i^X : 



(p) 
P o ■■ 



Xp) 



m« j_ 

A 4 + £ (f) 



(4 



(/>h2 



(43) 



w 22 _ A4 -r g( p) . 

We assume mechanical traction-free boundary conditions at infinity. Utilizing the jump 
conditions in the first of (7) and (8), along with the traction-free boundary conditions we have 
that 



d 



(b) _ >) 



< 2 -"2 4l=4l=^ (44) 

From the latter together with the second of Eq. (43) the expressions for the pressure in the 
phases are 

»/» dl 

(45) 



Pq — —a -r 



X L 



00' 



Substitution of Eq. (45) into the first of Eq. (43) leads to a relation between the applied La- 
grangian electric displacement field D = A 2 d and the resultant stretch ratio, namely 



l+Di 



1/6 



(46) 



where £ = yM/gW + v^/e^ an( * A = H^v^ + ji^v^. The latter appropriately reduces 
to the corresponding relation for a homogeneous specimen. In terms of the dimensionless 
quantities a = [i^/^ b \ j8 = /e^ b \ and D = D%/ y/jjffiefi), denoting the shear contrast, 
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permittivity contrast, and normalized Lagrangian electric displacement field, respectively, the 
stretch ratio is 

* = l+£ 2 ,„?_. ■ (47) 



v 



For the above choice of constitutive relation, we also have that Eq. (40) reduces to 

1 





= 4t = 


*g = 










- qg(p) - 


(?fP) - 


^121 


— ^211 _ 


^323 ~~ ' 


^1 111 


Op(/>) 

~ ®2121 — 


®2323 — 


^ 1212 


~ ^2222 ~~ 


®3232 — 



i(py 

w(p) _ 1 CSf(j>) _ 1 Ap) 



(48) 

d-!", (49) 



233 2 222 eip) 

<3 = <=f ( ^ 2 . (50) 

M (p) A 4 +^y4 P)2 - (51) 

The response of the laminate to a harmonic excitation superimposed on the aforementioned 
finite deformation is addressed next. Let x- and (p^ denote, respectively, the components 
of the incremental displacement and the incremental electric potential in the phases, such that 
e( p ) = — V<p( p \ The fields sought are to satisfy the incremental equations of motion along with 
Gauss equation (17) and the incompressibility constraint (18). Following Tiersten (1963), we 
consider fields that are functions of %2 and time alone, independent of x\ and X3, compatible 
with the simple-thickness mode (Mindlin, 1955). Consequently, in each phase the governing 
equations (17) and (18) simplify to 

A W *g2 = P W *g, (52a) 
^-^4 + P^)f^-P^-2^9^=P^ (52b) 



A W ^2 = P^S' (52c) 



and 



4 P 2 = 0, (53a) 
h$ = 0, (53b) 



respectively, where = /x^/X 4 , d { 2 p) = -e^ (d 2 <P ( 2 p) + ^drfi$\ Substitution of Eq. (53b) 



into Eq. (53a) reveals that <p^ = 0, hence the incremental electric potential is linear in xi, say 
<p( p ) = Lj +j^2 \ where l[ p ^ and are integration constants. More significantly, it 
implies that the incremental electrical and mechanical fields are not coupled, as Eq. (52b) re- 
mains a function of alone, independent of (p^ p \ together with Eq. (53a) which remains a 
function of <p^ p > alone, independent of h^ p \ Eq. (53b) states that along the thickness the incre- 
mental displacement is spatially constant. Together with Eqs. (52a) and (52c) the solution for 
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the incremental problem is 

x[ p) = (u { l p) smk^x 2 + U^ p) cosk^x 2 )e- 



iCOt 



xf = u( p) e- im , (54) 



= (u^hmk^x 2 + U [p) co%k^x 2 y- ia) \ 



where ki p ) = (0/ , = jl iP) / p ^ , CO is the angular excitation frequency, and Ur (r = 1,2,.. 
are integration constants. Without loss of generality the constants U-f can be set to be zero. 
The motion described by the solution Eq. (54) corresponds to displacement gradients which 
retain only shear components. For completeness, the solution for the incremental pressure 
emerges from Eq. (52b), which reduces to — 2 = P^ P ^2tv anc ^ i m pli es that pQ = pi p )e~' m , 
where pi p ) are integration constants. The incremental state vector of the n th cell takes the form 

e (p)_/*(p) M yiP) yiP) yiP) Jp) J(p)\t 
(«) ~~ I nl ' " 3 ' «12' Z 'n22' z 'n32' ^ ' w rc2 J • ^>~V 

The components of transfer matrix T^j relating the top sf^ and bottom sjjl state vectors are 
given in the Appendix. Note that the continuity conditions for x$ is identically satisfied when 
we set to zero. The corresponding solution of the eigenvalue problem Eq. (34) resulting 
from the Bloch-Floquet theorem can be expressed in terms of the trigonometric expression 

k B = \ arccos r\ , (56) 
h 



where 

1] = k 



cos^W fo W^ - \ (*«/*(-) +*<»>/*«) rin*W*>ain*W*« 



(57) 



with h^) = Hi p ) /A 2 and h = hi a ) + hi b ) . Eq. (56) is an extension of the solution given in Wang 
et al. (1986) to the class of finitely deformed infinite periodic DE laminates. The dependency 
on the bias field becomes evident when 77 is rephrased in terms of the electric bias field, the 
referential geometry and the frequency, namely 



' 2H V p. 



1/3 



aopi b )Hi b ) copi a )Hi a ) 



cos —7- cos 



api b ) pi°) \ copi b )Hi b ) copi a )Hi a 
+ -^-ttt sin — sin r 



p ui] : a P ih) r ^)(i + ^|) 1/3 M w(i+Dp^ 1/3 



(58) 



'2a 

In virtue of the uncoupled nature of the incremental problem, the incremental electric quanti- 
ties indeed do not appear in Eq. (58). We conclude this section noting that waves of certain 
frequencies cannot propagate when the pertained solution of ks is complex. These cases will 
be examined in the following section. 
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FIGURE 3. Variations of the Bloch wavenumber kg as functions of the nor- 
malized frequency a> = coH/c, for (a) a = /3 = 5, (b) a = /3 = 10 and (c) 
a = P = 20. The red, blue and black curves correspond to = 0.2, 0.5 and 
0.8, respectively. 
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FIGURE 4. Prohibited frequencies as functions of at D = 3, for (a) a = 
j8 = 5 and (b) a = j8 = 10 and (c) a = j8 = 20. 

4. Numerical examples 

With the help of a few examples we explore the dispersion relation characterizing the dy- 
namic response of the composite. We examine its dependency on the phases volume fractions 
and the contrast between the phases shear and dielectric moduli. In particular we are interested 
in the influence of the electric bias field on the dynamic response of the composite. To this end 
we make the following choice for the properties of phase b 



1000 



kg 



nr 



=200 [KPa], e, 



(b) 



(59) 



and determine dispersion diagrams for various choices of phase a properties and volume frac- 
tions at different values of the dimensionless bias field D. In all the forthcoming examples 



.(a) 



p^ is assumed. 



Fig. 3 shows variations of the Bloch wavenumber kg as functions of the normalized fre- 
quency G) = coH/c where c = y/jl/p, p = p^v" + p^v", for a representative value of the 
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Figure 5. Prohibited frequencies as functions of a at D = 3 and = 0.5, 
for (a) j8 =0.1 and (b) j8 = 10. 

bias field D = 3. Figs. 3a, b and c correspond toct = /3 = 5,10 and 20, respectively. The red, 
blue and black curves correspond to v^) = 0.2,0.5 and 0.8, respectively. The period of the 
dispersion curves becomes smaller with an increase of the softer phase volume fraction. An 
inverse effect is revealed when the contrast is enhanced, that is longer periods for higher values 
of a and /3. 

Fig. 4 displays the regions of prohibited normalized frequencies G) as functions of the matrix 
volume fraction at a fixed value of the bias field D = 3. Figs. 4a, b and c correspond to 
a = /3 = 5, 10 and 20, respectively. Herein and henceforth the band-gaps are the blue regions in 
the plots. We observe how an increase of the softer phase volume fraction results in appearance 
of additional thinner bands. Conversely, when the contrasts are increased, that is higher values 
of a and /3, there are fewer thicker bands. The latter is evident in Fig. 5 which displays the 
prohibited frequencies as functions of a for = 0.5 and D = 3, with different values of 
/3. Specifically, Fig. 5a corresponds to /3 = 0.1 and Fig. 5b corresponds to j8 = 10. Similar 
trends are observed when the dielectric contrast is investigated. Fig. 6 displays the prohibited 
frequencies (b as functions of /3 for v'*' = 0.5, with different values of a. Fig. 6a corresponds 
to a = 0.1 and Fig. 6b to a = 10. Figs. 5 and 6 reveal that whenever the phases have non- 
monotonous ratios of shear moduli and dielectric contrasts, i.e. a > 1 and /3 > 1 or a < 1 and 
/3 > 1, the band gaps become thiner and denser. 

Fig. 7 illustrates the influence of the bias field by examining the variations of the BGs as 
functions of the normalized parameter D. Values of v^ b ' =0.8 and a = /3 = 5 were chosen in 
Fig. 7a, = 0.5 and a = j8 = 10 in Fig. 7b, and = 0.2 and a = j8 = 20 in Fig. 7c. The 
influence of the electric field is evident as the BGs are shifted toward lower frequencies when 
the field is increased. This phenomena suggests that control of BGs and filtering of desired 
frequencies are feasible by activating DE laminates with adequate bias electric fields. 
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FIGURE 6. Prohibited frequencies as functions of j8 at D = 3 and = 0.5, 
for (a) a = 0.1 and (b) a = 10. 
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Figure 7. Prohibited frequencies as functions of D for (a) = 0.8, a = j8 = 
5, (b) vW = 0.5, a = j8 = 10 and (c) = 0.2, a = j8 = 20. 

Finally, we consider realizable laminates made out of commercially available DEs. We ex- 
amine three materials: VHB-4910 by 3M, ELASTOSIL RT-625 by Wacker, and fluorosilicone 
730 by Dow Corning. Approximate physical properties of these elastomers, as reported in the 
literature (e.g., Kofod and Sommer-Larsen, 2005; Kornbluh and Pelrine, 2008), are summa- 
rized in table 1. We investigate different combinations of the two-phase laminate and assume 
thatv^ =vH 

Fig. 8 displays the dispersion relation in terms of Bloch wavenumber kg as a function of the 
normalized frequency 6). Fig. 8a corresponds to a combination of VHB-4910 and ELASTOSIL 



Table 1 . Physical properties of commercially available DEs. 



Material 


density [kg/m 3 ] 


shear modulus [kPa] 


relative dielectric 
constant 


VHB-4910 


960 


406 


4.7 


ELASTOSIL RT-625 


1020 


342 


2.7 


Fluorosilicone 730 


1400 


167 


6.9 
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VHB-4910 and ELASTOSIL RT-625 VHB-4910 and fluorosilicone 730 




Figure 8. Variations of the Bloch wavenumber kg as functions of the nor- 
malized frequency 6) = CoH jc for combinations of VHB-4910 with (a) ELAS- 
TOSIL RT-625, and with (b) fluorosilicone 730. 
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Figure 9. Prohibited frequencies as functions of D for combinations of VHB- 
4910 with (a) ELASTOSIL RT-625, and (b) fluorosilicone 730. 

RT-625, while Fig. 8b to VHB-4910 and fluorosilicone 730. In terms of the dimensionless 
parameters in Fig. 8a a = 1.19 and /3 = 1 .74, and in Fig. 8b a = 2.44 and /3 = 0.68. The blue, 
black and red curves correspond to D = 0.5, 1 and 1.5, respectively. 

Fig. 9 shows the prohibited frequencies as functions of the bias field in terms of the di- 
mensionless bias field D. Again, Fig. 9a corresponds to the combination of VHB-4910 and 
ELASTOSIL RT-625, and Fig.9b c to VHB-4910 and fluorosilicone 730. Wider bands are de- 
picted in Fig. 9b in comparison with Fig. 9a, since the latter corresponds to a laminate with 
higher shear contrast between the phases. 

5. Concluding remarks 

Motivated by the ability of DEs to undergo large deformations and change their mechanical 
and electrical properties when electrostatically excited, the feasibility of inducing and con- 
trolling BGs by the electrical bias field is explored. To this end we first considered the static 
response of an infinite periodic laminate to a bias electric displacement field. Subsequently, 
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the propagation of small waves superposed on the pre-strained composite is addressed. Appli- 
cations of the transfer matrix method along with the Bloch-Floquet theorem resulted in a dis- 
persion relation between the waves frequencies, the phase velocities, and the Bloch-parameter. 
Numerical analysis of these relations was conducted for specific two-phase laminates in order 
to examine how the stop-bands, associated with complex Bloch-parameter, vary as functions 
of the geometrical, mechanical and electric properties of the laminate. Most importantly, this 
analysis reveals the influence of the electrostatic bias fields on the BGs. Our findings demon- 
strate how when the concentration of the softer phase in increased, the bands become thinner 
and additional band gaps the appear. An inverse effect is observed when the contrast between 
the phases properties is increased (larger values of a and /3). The primary conclusion of our 
analysis is depicted in Fig. 7, demonstrating how variations in the bias electric displacement 
field lead to modifications in the width and shifts in the range of the prohibited frequencies. 
We note that for lower values of a and /3 the effect of the bias electric field becomes more 
evident in the sense that the shift in the BGs range is larger. From practical viewpoint point, 
while stopping larger bands of frequencies is feasible by choosing a higher contrast, lower 
contrasts will allow to actively adjust the desired band at a higher precision. The analysis was 
concluded with two examples in which we chose phases properties to be identical to those 
of commercially available DEs. The results show how tunable stop-bands are achievable by 
properly adjusting the electrostatic bias field. 

Appendix A. The transfer matrix 
The components of the transfer matrices T^j appearing in section 3: 

r/f > = T$ = 4 p) = 4 P) = cos kWhW , (60) 

= t£ ] = sinJb^kO /fik^\ (61) 

T 3 ( f } = T$ = T 5 f = -p.kW sin^W^, (62) 

T$ = -h^/de (p \ (63) 

T$ = T$ = l, (64) 

with all the remaining components being zero. 
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